Symbiosis: In search of a deeper understanding.

How do distinct species cofunction in symbiosis, despite conflicting interests? A new collection of articles explores emerging themes as researchers exploit modern research tools and new models to unravel how symbiotic interactions function and evolve.

Duplication and neofunctionalization of a horizontally-transferred xyloglucanase as a facet of the red queen co-evolutionary dynamic.

Oomycetes are heterotrophic protists that share phenotypic similarities with fungi, including the ability to cause plant diseases, but branch in a separate and distant region of the eukaryotic tree of life. It has been suggested that multiple horizontal gene transfers (HGTs) from fungi-to-oomycetes contributed to the evolution of plant-pathogenic traits. These HGTs are predicted to include secreted proteins that degrade plant cell walls. This is a key trait in the pathology of many oomycetes, as the plant cell wall represents a primary barrier to pathogen invasion and a rich source of carbohydrates. Many of the HGT gene families identified have undergone multiple rounds of duplication. Using a combination of phylogenomic analysis and functional assays, we investigate the diversification of a horizontally-transferred xyloglucanase gene family in the model oomycete species Phytophthora sojae. Our analyses detect 11 genes retained in P. sojae among a complex pattern of gene duplications and losses. Using a phenotype assay, based on heterologous expression in yeast, we show that eight of these paralogs have xyloglucanase function, including variants with distinct protein characteristics, such as a long-disordered C-terminal extension that can increase xyloglucanase activity. The functional xyloglucanase variants analysed subtend an ancestral node close to the fungi-oomycetes gene transfer, suggesting the horizontally-transferred gene was a bona fide xyloglucanase. Expression of xyloglucanase paralogs in Nicotiana benthamiana triggers distinct patterns of reactive oxygen species (ROS) generation, demonstrating that enzyme variants differentially stimulate pattern-triggered immunity in plants. Mass spectrometry of detectable enzymatic products demonstrates that some paralogs catalyze production of variant breakdown profiles, suggesting that secretion of multiple xyloglucanase variants increases efficiency of xyloglucan breakdown, as well as potentially diversifying the range of Damage-Associated Molecular Patterns (DAMPs) released during pathogen attack. We suggest that such patterns of protein neofunctionalization, and variant host responses, represent an aspect of the Red Queen host-pathogen co-evolutionary dynamic. Significance Statement: The oomycetes are a diverse group of eukaryotic microbes that include some of the most devastating pathogens of plants. Oomycetes perceive, invade, and colonize plants in similar ways to fungi, in part because they acquired the genes to attack and feed on plants from fungi. These genes are predicted to be useful to oomycete plant pathogens because they have undergone multiple rounds of gene duplication. One key enzyme for attacking plant cell wall structures is called xyloglucanase. Xyloglucanase in the oomycetes has undergone multiple rounds of gene duplication, leading to variants including an enzyme with a C-terminal extension that increases activity. Some xyloglucanase variants trigger unique patterns of reactive oxygen species (ROS) in planta, and generate different profiles of cell wall breakdown products - such outcomes could act to mystify and increase the workload of the plant immune system, allowing successful pathogens to proliferate.

Identification of a non-canonical ciliate nuclear genetic code where UAA and UAG code for different amino acids.

The genetic code is one of the most highly conserved features across life. Only a few lineages have deviated from the "universal" genetic code. Amongst the few variants of the genetic code reported to date, the codons UAA and UAG virtually always have the same translation, suggesting that their evolution is coupled. Here, we report the genome and transcriptome sequencing of a novel uncultured ciliate, belonging to the Oligohymenophorea class, where the translation of the UAA and UAG stop codons have changed to specify different amino acids. Genomic and transcriptomic analyses revealed that UAA has been reassigned to encode lysine, while UAG has been reassigned to encode glutamic acid. We identified multiple suppressor tRNA genes with anticodons complementary to the reassigned codons. We show that the retained UGA stop codon is enriched in the 3'UTR immediately downstream of the coding region of genes, suggesting that there is functional drive to maintain tandem stop codons. Using a phylogenomics approach, we reconstructed the ciliate phylogeny and mapped genetic code changes, highlighting the remarkable number of independent genetic code changes within the Ciliophora group of protists. According to our knowledge, this is the first report of a genetic code variant where UAA and UAG encode different amino acids.

Thomas A. Richards.

Interview with Tom Richards, who uses comparative genomics and molecular biology to explore eukaryotic cell evolution.

Systematic comparison of unilamellar vesicles reveals that archaeal core lipid membranes are more permeable than bacterial membranes.

One of the deepest branches in the tree of life separates the Archaea from the Bacteria. These prokaryotic groups have distinct cellular systems including fundamentally different phospholipid membrane bilayers. This dichotomy has been termed the lipid divide and possibly bestows different biophysical and biochemical characteristics on each cell type. Classic experiments suggest that bacterial membranes (formed from lipids extracted from Escherichia coli, for example) show permeability to key metabolites comparable to archaeal membranes (formed from lipids extracted from Halobacterium salinarum), yet systematic analyses based on direct measurements of membrane permeability are absent. Here, we develop a new approach for assessing the membrane permeability of approximately 10 µm unilamellar vesicles, consisting of an aqueous medium enclosed by a single lipid bilayer. Comparing the permeability of 18 metabolites demonstrates that diether glycerol-1-phosphate lipids with methyl branches, often the most abundant membrane lipids of sampled archaea, are permeable to a wide range of compounds useful for core metabolic networks, including amino acids, sugars, and nucleobases. Permeability is significantly lower in diester glycerol-3-phosphate lipids without methyl branches, the common building block of bacterial membranes. To identify the membrane characteristics that determine permeability, we use this experimental platform to test a variety of lipid forms bearing a diversity of intermediate characteristics. We found that increased membrane permeability is dependent on both the methyl branches on the lipid tails and the ether bond between the tails and the head group, both of which are present on the archaeal phospholipids. These permeability differences must have had profound effects on the cell physiology and proteome evolution of early prokaryotic forms. To explore this further, we compare the abundance and distribution of transmembrane transporter-encoding protein families present on genomes sampled from across the prokaryotic tree of life. These data demonstrate that archaea tend to have a reduced repertoire of transporter gene families, consistent with increased membrane permeation. These results demonstrate that the lipid divide demarcates a clear difference in permeability function with implications for understanding some of the earliest transitions in cell origins and evolution.

Transporter Proteins as Ecological Assets and Features of Microbial Eukaryotic Pangenomes.

Here we review two connected themes in evolutionary microbiology: (a) the nature of gene repertoire variation within species groups (pangenomes) and (b) the concept of metabolite transporters as accessory proteins capable of providing niche-defining "bolt-on" phenotypes. We discuss the need for improved sampling and understanding of pangenome variation in eukaryotic microbes. We then review the factors that shape the repertoire of accessory genes within pangenomes. As part of this discussion, we outline how gene duplication is a key factor in both eukaryotic pangenome variation and transporter gene family evolution. We go on to outline how, through functional characterization of transporter-encoding genes, in combination with analyses of how transporter genes are gained and lost from accessory genomes, we can reveal much about the niche range, the ecology, and the evolution of virulence of microbes. We advocate for the coordinated systematic study of eukaryotic pangenomes through genome sequencing and the functional analysis of genes found within the accessory gene repertoire.

Shuffled ATG8 interacting motifs form an ancestral bridge between UFMylation and autophagy.

UFMylation involves the covalent modification of substrate proteins with UFM1 (Ubiquitin-fold modifier 1) and is important for maintaining ER homeostasis. Stalled translation triggers the UFMylation of ER-bound ribosomes and activates C53-mediated autophagy to clear toxic polypeptides. C53 contains noncanonical shuffled ATG8-interacting motifs (sAIMs) that are essential for ATG8 interaction and autophagy initiation. However, the mechanistic basis of sAIM-mediated ATG8 interaction remains unknown. Here, we show that C53 and sAIMs are conserved across eukaryotes but secondarily lost in fungi and various algal lineages. Biochemical assays showed that the unicellular alga Chlamydomonas reinhardtii has a functional UFMylation pathway, refuting the assumption that UFMylation is linked to multicellularity. Comparative structural analyses revealed that both UFM1 and ATG8 bind sAIMs in C53, but in a distinct way. Conversion of sAIMs into canonical AIMs impaired binding of C53 to UFM1, while strengthening ATG8 binding. Increased ATG8 binding led to the autoactivation of the C53 pathway and sensitization of Arabidopsis thaliana to ER stress. Altogether, our findings reveal an ancestral role of sAIMs in UFMylation-dependent fine-tuning of C53-mediated autophagy activation.

Amphibian Perkinsea.

Vanessa Smilansky and Thomas A. Richards introduce Perkinsea - a lineage of freshwater parasitic protists that infect certain amphibians and cause of severe Perkinsea infection.

A Genome Sequence Assembly of the Phototactic and Optogenetic Model Fungus Blastocladiella emersonii Reveals a Diversified Nucleotide-Cyclase Repertoire.

The chytrid fungus Blastocladiella emersonii produces spores with swimming tails (zoospores); these cells can sense and swim toward light. Interest in this species stems from ongoing efforts to develop B. emersonii as a model for understanding the evolution of phototaxis and the molecular cell biology of the associated optogenetic circuits. Here, we report a highly contiguous genome assembly and gene annotation of the B. emersonii American Type Culture Collection 22665 strain. We integrate a PacBio long-read library with an Illumina paired-end genomic sequence survey leading to an assembly of 21 contigs totaling 34.27 Mb. Using these data, we assess the diversity of sensory system encoding genes. These analyses identify a rich complement of G-protein-coupled receptors, ion transporters, and nucleotide cyclases, all of which have been diversified by domain recombination and tandem duplication. In many cases, these domain combinations have led to the fusion of a protein domain to a transmembrane domain, tying a putative signaling function to the cell membrane. This pattern is consistent with the diversification of the B. emersonii sensory-signaling systems, which likely plays a varied role in the complex life cycle of this fungus.

A diversified and segregated mRNA spliced-leader system in the parasitic Perkinsozoa.

Spliced-leader trans-splicing (SLTS) has been described in distantly related eukaryotes and acts to mark mRNAs with a short 5' exon, giving different mRNAs identical 5' sequence-signatures. The function of these systems is obscure. Perkinsozoa encompasses a diversity of parasitic protists that infect bivalves, toxic-tide dinoflagellates, fish and frog tadpoles. Here, we report considerable sequence variation in the SLTS-system across the Perkinsozoa and find that multiple variant SLTS-systems are encoded in parallel in the ecologically important Perkinsozoa parasite Parvilucifera sinerae. These results demonstrate that the transcriptome of P. sinerae is segregated based on the addition of different spliced-leader (SL) exons. This segregation marks different gene categories, suggesting that SL-segregation relates to functional differentiation of the transcriptome. By contrast, both sets of gene categories are present in the single SL-transcript type sampled from Maranthos, implying that the SL-segregation of the Parvilucifera transcriptome is a recent evolutionary innovation. Furthermore, we show that the SLTS-system marks a subsection of the transcriptome with increased mRNA abundance and includes genes that encode the spliceosome system necessary for SLTS-function. Collectively, these data provide a picture of how the SLTS-systems can vary within a major evolutionary group and identify how additional transcriptional-complexity can be achieved through SL-segregation.

Exaggerated trans-membrane charge of ammonium transporters in nutrient-poor marine environments.

Transporter proteins are a vital interface between cells and their environment. In nutrient-limited environments, microbes with transporters that are effective at bringing substrates into their cells will gain a competitive advantage over variants with reduced transport function. Microbial ammonium transporters (Amt) bring ammonium into the cytoplasm from the surrounding periplasm space, but diagnosing Amt adaptations to low nutrient environments solely from sequence data has been elusive. Here, we report altered Amt sequence amino acid distribution from deep marine samples compared to variants sampled from shallow water in two important microbial lineages of the marine water column community-Marine Group I Archaea (Thermoproteota) and the uncultivated gammaproteobacterial lineage SAR86. This pattern indicates an evolutionary pressure towards an increasing dipole in Amt for these clades in deep ocean environments and is predicted to generate stronger electric fields facilitating ammonium acquisition. This pattern of increasing dipole charge with depth was not observed in lineages capable of accessing alternative nitrogen sources, including the abundant alphaproteobacterial clade SAR11. We speculate that competition for ammonium in the deep ocean drives transporter sequence evolution. The low concentration of ammonium in the deep ocean is therefore likely due to rapid uptake by Amts concurrent with decreasing nutrient flux.

A light-sensing system in the common ancestor of the fungi.

Diverse light-sensing organs (i.e., eyes) have evolved across animals. Interestingly, several subcellular analogs have been found in eukaryotic microbes.1 All of these systems have a common "recipe": a light occluding or refractory surface juxtaposed to a membrane-layer enriched in type I rhodopsins.1-4 In the fungi, several lineages have been shown to detect light using a diversity of non-homologous photo-responsive proteins.5-7 However, these systems are not associated with an eyespot-like organelle with one exception found in the zoosporic fungus Blastocladiella emersonii (Be).8Be possesses both elements of this recipe: an eyespot composed of lipid-filled structures (often called the side-body complex [SBC]), co-localized with a membrane enriched with a gene-fusion protein composed of a type I (microbial) rhodopsin and guanylyl cyclase enzyme domain (CyclOp-fusion protein).8,9 Here, we identify homologous pathway components in four Chytridiomycota orders (Chytridiales, Synchytriales, Rhizophydiales, and Monoblepharidiales). To further explore the architecture of the fungal zoospore and its lipid organelles, we reviewed electron microscopy data (e.g., the works of Barr and Hartmann10 and Reichle and Fuller11) and performed fluorescence-microscopy imaging of four CyclOp-carrying zoosporic fungal species, showing the presence of a variety of candidate eyespot-cytoskeletal ultrastructure systems. We then assessed the presence of canonical photoreceptors across the fungi and inferred that the last common fungal ancestor was able to sense light across a range of wavelengths using a variety of systems, including blue-green-light detection. Our data imply, independently of how the fungal tree of life is rooted, that the apparatus for a CyclOp-organelle light perception system was an ancestral feature of the fungi.

Nutrient and salt depletion synergistically boosts glucose metabolism in individual Escherichia coli cells.

The interaction between a cell and its environment shapes fundamental intracellular processes such as cellular metabolism. In most cases growth rate is treated as a proximal metric for understanding the cellular metabolic status. However, changes in growth rate might not reflect metabolic variations in individuals responding to environmental fluctuations. Here we use single-cell microfluidics-microscopy combined with transcriptomics, proteomics and mathematical modelling to quantify the accumulation of glucose within Escherichia coli cells. In contrast to the current consensus, we reveal that environmental conditions which are comparatively unfavourable for growth, where both nutrients and salinity are depleted, increase glucose accumulation rates in individual bacteria and population subsets. We find that these changes in metabolic function are underpinned by variations at the translational and posttranslational level but not at the transcriptional level and are not dictated by changes in cell size. The metabolic response-characteristics identified greatly advance our fundamental understanding of the interactions between bacteria and their environment and have important ramifications when investigating cellular processes where salinity plays an important role.

Systematic evaluation of horizontal gene transfer between eukaryotes and viruses.

Gene exchange between viruses and their hosts acts as a key facilitator of horizontal gene transfer and is hypothesized to be a major driver of evolutionary change. Our understanding of this process comes primarily from bacteria and phage co-evolution, but the mode and functional importance of gene transfers between eukaryotes and their viruses remain anecdotal. Here we systematically characterized viral-eukaryotic gene exchange across eukaryotic and viral diversity, identifying thousands of transfers and revealing their frequency, taxonomic distribution and projected functions. Eukaryote-derived viral genes, abundant in the Nucleocytoviricota, highlighted common strategies for viral host-manipulation, including metabolic reprogramming, proteolytic degradation and extracellular modification. Furthermore, viral-derived eukaryotic genes implicate genetic exchange in the early evolution and diversification of eukaryotes, particularly through viral-derived glycosyltransferases, which have impacted structures as diverse as algal cell walls, trypanosome mitochondria and animal tissues. These findings illuminate the nature of viral-eukaryotic gene exchange and its impact on the evolution of viruses and their eukaryotic hosts.

A cell-cell atlas approach for understanding symbiotic interactions between microbes.

Natural environments are composed of a huge diversity of microorganisms interacting with each other to form complex functional networks. Our understanding of the operative nature of host-symbiont associations is limited because propagating such associations in a laboratory is challenging. The advent of single-cell technologies applied to, for example, animal cells and apicomplexan parasites has revolutionized our understanding of development and disease. Such cell atlas approaches generate maps of cell-specific processes and variations within cellular populations. These methods can now be combined with cellular-imaging so that interaction stage versus transcriptome state can be quantized for microbe-microbe interactions. We predict that the combination of these methods applied to the study of symbioses will transform our understanding of many ecological interactions, including those sampled directly from natural environments.

Emergent RNA-RNA interactions can promote stability in a facultative phototrophic endosymbiosis.

Eukaryote-eukaryote endosymbiosis was responsible for the spread of chloroplast (plastid) organelles. Stability is required for the metabolic and genetic integration that drives the establishment of new organelles, yet the mechanisms that act to stabilize emergent endosymbioses-between two fundamentally selfish biological organisms-are unclear. Theory suggests that enforcement mechanisms, which punish misbehavior, may act to stabilize such interactions by resolving conflict. However, how such mechanisms can emerge in a facultative endosymbiosis has yet to be explored. Here, we propose that endosymbiont-host RNA-RNA interactions, arising from digestion of the endosymbiont population, can result in a cost to host growth for breakdown of the endosymbiosis. Using the model facultative endosymbiosis between Paramecium bursaria and Chlorella spp., we demonstrate that this mechanism is dependent on the host RNA-interference (RNAi) system. We reveal through small RNA (sRNA) sequencing that endosymbiont-derived messenger RNA (mRNA) released upon endosymbiont digestion can be processed by the host RNAi system into 23-nt sRNA. We predict multiple regions of shared sequence identity between endosymbiont and host mRNA, and demonstrate through delivery of synthetic endosymbiont sRNA that exposure to these regions can knock down expression of complementary host genes, resulting in a cost to host growth. This process of host gene knockdown in response to endosymbiont-derived RNA processing by host RNAi factors, which we term "RNAi collisions," represents a mechanism that can promote stability in a facultative eukaryote-eukaryote endosymbiosis. Specifically, by imposing a cost for breakdown of the endosymbiosis, endosymbiont-host RNA-RNA interactions may drive maintenance of the symbiosis across fluctuating ecological conditions.

Expanded host and geographic range of tadpole associations with the Severe Perkinsea Infection group.

Severe Perkinsea infection is an emerging disease of amphibians, specifically tadpoles. Disease presentation correlates with liver infections of a subclade of Perkinsea (Alveolata) protists, named Pathogenic Perkinsea Clade (PPC). Tadpole mortality events associated with PPC infections have been reported across North America, from Alaska to Florida. Here, we investigate the geographic and host range of PPC associations in seemingly healthy tadpoles sampled from Panama, a biogeographic provenance critically affected by amphibian decline. To complement this work, we also investigate a mortality event among Hyla arborea tadpoles in captive-bred UK specimens. PPC SSU rDNA was detected in 10 of 81 Panama tadpoles tested, and H. arborea tadpoles from the UK. Phylogenies of the Perkinsea SSU rDNA sequences demonstrate they are highly similar to PPC sequences sampled from mortality events in the USA, and phylogenetic analysis of tadpole mitochondrial SSU rDNA demonstrates, for the first time, PPC associations in diverse hylids. These data provide further understanding of the biogeography and host range of this putative pathogenic group, factors likely to be important for conservation planning.

Characterization of the RNA-interference pathway as a tool for reverse genetic analysis in the nascent phototrophic endosymbiosis, Paramecium bursaria.

Endosymbiosis was fundamental for the evolution of eukaryotic complexity. Endosymbiotic interactions can be dissected through forward- and reverse-genetic experiments, such as RNA-interference (RNAi). However, distinguishing small (s)RNA pathways in a eukaryote-eukaryote endosymbiotic interaction is challenging. Here, we investigate the repertoire of RNAi pathway protein-encoding genes in the model nascent endosymbiotic system, Paramecium bursaria-Chlorella spp. Using comparative genomics and transcriptomics supported by phylogenetics, we identify essential proteome components of the small interfering (si)RNA, scan (scn)RNA and internal eliminated sequence (ies)RNA pathways. Our analyses reveal that copies of these components have been retained throughout successive whole genome duplication (WGD) events in the Paramecium clade. We validate feeding-induced siRNA-based RNAi in P. bursaria via knock-down of the splicing factor, u2af1, which we show to be crucial to host growth. Finally, using simultaneous knock-down 'paradox' controls to rescue the effect of u2af1 knock-down, we demonstrate that feeding-induced RNAi in P. bursaria is dependent upon a core pathway of host-encoded Dcr1, Piwi and Pds1 components. Our experiments confirm the presence of a functional, host-derived RNAi pathway in P. bursaria that generates 23-nt siRNA, validating the use of the P. bursaria-Chlorella spp. system to investigate the genetic basis of a nascent endosymbiosis.

Single-cell genomics unveils a canonical origin of the diverse mitochondrial genomes of euglenozoans.

BACKGROUND: The supergroup Euglenozoa unites heterotrophic flagellates from three major clades, kinetoplastids, diplonemids, and euglenids, each of which exhibits extremely divergent mitochondrial characteristics. Mitochondrial genomes (mtDNAs) of euglenids comprise multiple linear chromosomes carrying single genes, whereas mitochondrial chromosomes are circular non-catenated in diplonemids, but circular and catenated in kinetoplastids. In diplonemids and kinetoplastids, mitochondrial mRNAs require extensive and diverse editing and/or trans-splicing to produce mature transcripts. All known euglenozoan mtDNAs exhibit extremely short mitochondrial small (rns) and large (rnl) subunit rRNA genes, and absence of tRNA genes. How these features evolved from an ancestral bacteria-like circular mitochondrial genome remains unanswered. RESULTS: We sequenced and assembled 20 euglenozoan single-cell amplified genomes (SAGs). In our phylogenetic and phylogenomic analyses, three SAGs were placed within kinetoplastids, 14 within diplonemids, one (EU2) within euglenids, and two SAGs with nearly identical small subunit rRNA gene (18S) sequences (EU17/18) branched as either a basal lineage of euglenids, or as a sister to all euglenozoans. Near-complete mitochondrial genomes were identified in EU2 and EU17/18. Surprisingly, both EU2 and EU17/18 mitochondrial contigs contained multiple genes and one tRNA gene. Furthermore, EU17/18 mtDNA possessed several features unique among euglenozoans including full-length rns and rnl genes, six mitoribosomal genes, and nad11, all likely on a single chromosome. CONCLUSIONS: Our data strongly suggest that EU17/18 is an early-branching euglenozoan with numerous ancestral mitochondrial features. Collectively these data contribute to untangling the early evolution of euglenozoan mitochondria.

A novel duplex qPCR assay for stepwise detection of multiple Perkinsea protistan infections of amphibian tissues.

Alveolate protists within the phylum Perkinsea have been found to infect amphibians across a broad taxonomic and geographic range. Phylogenetic analysis has suggested the existence of two clades of amphibian Perkinsea: a putatively non-pathogenic clade linked to asymptomatic infections of tadpoles in Africa, Europe and South America, and a putatively pathogenic clade linked to disease and mass mortality events of tadpoles in North America. Here, we describe the development of a duplex TaqMan qPCR assay to detect and discriminate between rDNA sequences from both clades of Perkinsea in amphibian tissues. The assay uses a single primer pair to target an 18S small subunit (SSU) ribosomal RNA (rRNA) gene region shared between the two clades, and two dual-labelled probes to target a region within this fragment that is diagnostic for each clade. This assay enables rapid screening for each of the two Perkinsea groups, allowing for detection, primarily of the phylogenetic group associated with disease outbreaks, and secondarily for the phylogenetic group with no current disease relationship identified. Incorporation of our novel qPCR assay into the routine surveillance of amphibian populations will allow for the assessment of the incidence of each protist clade, thereby providing an improved understanding of Perkinsea infection pervasiveness and a method to underpin future conservation planning.